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We examine the spatial distribution of electrons generated by a fixed energy point source in 
uniform, parallel electric and magnetic fields. This problem is simple enough to permit analytic 
quantum and semiclassical solution, and it harbors a rich set of features which find their interpreta- 
tion in the unusual and interesting properties of the classical motion of the electrons; For instance, 
the number of interfering trajectories can be adjusted in this system, and the turning surfaces of 
classical motion contain a complex array of singularities. We perform a comprehensive analysis 
of both the semiclassical approximation and the quantum solution, and we make predictions that 
should serve as a guide for future photodetachment experiments. 

PACS numbers: 03.65.Sq, 03.75.-b, 32.80.Gc 



I. INTRODUCTION 

The photoabsorption cross-section of atoms or nega- 
tive ions in presence of external electric 1', '5, 'sj or mag- 
netic 1^] fields shows oscillations as a function of the 
photon energy. For example, in a recent experiment, Yu- 
kich et al. Q measured the photodetachment rates of 
negative ions in a magnetic field and showed that the 
addition of a parallel electric field causes variations in 
the cross-section. Closed Orbit Theory interprets the 
oscillations as interferences between waves that go out 
from, and return to, the atomic source 15 lol llolllll|. 
Such semiclassical models, as well as quantum calcula- 
tions 0, 0, 0, 0, 0| , are in good agreement with ex- 
perimental results. 

Many years ago, Demkov, Kondratovich and Ostro- 
vskii |l3 suggested that interference structures could also 
be observed in a spatially resolved measurement of the 
detached electrons: If electrons of fixed energy emerge 
from a point source, and these electrons are accelerated 
by a uniform electric field toward a microchannelplate 
detector, then two trajectories go from the source to 
any given point on the detector |17|. Along each tra- 
jectory the wave function accumulates a phase propor- 
tional to the classical action of that path. The re- 
sulting electron distribution at the detector shows a 
set of concentric interference fringes, associated with 
the two paths linking the source with the detection 
point. More complicated paths and associated inter- 
ference structures arise in photoionization experiments 
[T5 Esl I20I I21I . The observation and interpretation of 
such interference patterns is called photodetachment or 
photoionization microsco py. These techniques were pio- 
neered by Blondel et al. p^. l23l|. and Bordas, Vrakking 
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and co-worke rs 25l , respectively. Quantum cal- 
culations 0, Il6l I27 1 yield excellent agreement with the 
experimental data, and in turn allow the determination 
of electron affinities from the interference patterns with 
unprecedented accuracy . 

In this paper we predict the distribution of electrons 
that would be seen in an ideal photodetachment micro- 
scope placed in parallel, uniform, static electric and mag- 
netic fields 6 and B. With the idealization that the 
electrons come from a point source at fixed energy, we 
compute the spatial distribution of electrons, and pre- 
dict their current density arriving at various points on a 
detector (Fig.nj. 

A weak cw laser field acting on a negative ion gives 
rise to outgoing, coherent, fixed-energy electron waves 
(compare a vibrator in a ripple tank). The angular dis- 
tribution of these waves is governed by selection rules 
for the initial and final state; for simplici ty, w e assume a 
spherically symmetric emission pattern |53- The out- 
ward propagation of these waves can be described by 
considering the evolution of the corresponding radially 
outgoing classical trajectories. As the paths bend in the 
applied electric and magnetic fields, the waves refract to 
follow them. For weak external fields, the classical range 
of motion of the ejected electrons extends far beyond the 
residual atomic core, so their interaction with the core is 
negligible, and rescattering is limited to a small fraction 
of the initial outgoing wave. It is a good approximation 
to neglect the presence of the atom after the photodetach- 
ment event altogether. Then the dynamics of the elec- 
trons is governed solely by the externally applied fields. 
(Rescattering is addressed e. g. in Ref. 30., 31. .32j.) 

In parallel fields aligned with the z axis, every elec- 
tron undergoes uniformly accelerated motion along the 
symmetry axis while orbiting on a circular track in the 
perpendicular x — y plane (Section^. The initial emis- 
sion angle 0' of the orbit relative to the field direction 
fixes both the energy contained in the motion along the 
z axis and the energy available for transverse cyclotron 
motion, and these two energies are separately conserved. 
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FIG. 1: Motion of electrons emitted from a fixed energy, isotropic point source at the origin into parallel fields. (Parameters 
used: E = 10"* eV, £ = f5V/m, B = 0.02 T, corresponding to a force parameter 77 = voB/£ = 7.908 0. All dimensions 
are in fim.) The classical trajectory field (top left panel, blue curves) traces out an intersecting set of caustic surfaces (red 
curves) that separate regions with different numbers of classical paths connecting the source to a given destination point. The 
maximum number (here, eight) occurs in the diamond-shaped entities along the symmetry axis. The trajectory field contains 
two closed orbits that return to the source fio| . a "snake" orbit (green), and a "balloon" orbit (purple). — The center panel 
depicts the density profile of the electrons in primitive semiclassical approximation (blue plot), superimposed onto the caustic 
surfaces (red) and the closed orbits. The caustics clearly delineate the features of the density distribution, and the interference 
pattern shows approximate symmetry with respect to the closed "snake" orbit. (The letter combinations (green) refer to the 
classification of the caustics in Section UTTl ) — The quantum density profile is displayed in the top right panel. Its agreement 
with the semiclassical result (center panel) is remarkable. — Bottom row: Profiles of the radially integrated electronic density 
distribution n{p) at a distance z = 30 /im from the source (indicated by purple arrows above) for a normalized fiux of one 
electron per second, in units of m~^. The green curve in the left panel denotes the purely classical cross section. It is strongly 
modified by semiclassical interference (center panel). Both classical and semiclassical density profiles diverge at the positions 
of the caustics (marked by red arrows) . Otherwise, the semiclassical density profile quantitatively matches the quantum result 
(blue curve, right panel). 
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The classical spatial distribution is determined by the set 
of all trajectories, combined with the angular characteris- 
tics of the source (isotropic here). The upper left panel of 
Fig-Ddisplays a profile of the resulting trajectory pattern 
(blue) for E = 15V/m, B ^ 0.02 T, and E = lO-^eV. 
(The ordinate denotes the radial distance p from the z 
axis. All dimensions are in /im; axes are not to scale.) 
The trajectories trace out an infinite sequence of caus- 
tic surfaces (red curves); if we move a detector across 
a caustic surface, the number of classical paths arriving 
at the detector changes. These caustic surfaces usually 
begin and end on the symmetry axis in a rotationally 
symmetric cusp. We will return to the caustic surfaces 
and their singularities in Section Hill 

The lower left plot of Fig. ^ shows a profile of the in- 
tegrated classical electron density profile nci{p) = 2Trp ■ 
Pc\{p) (green), evaluated at a distance z — 30 /im from 
the source. It is marked by divergences at the inter- 
sections with the caustic surfaces. This classical result 
reproduces the overall trends in the quantum distribu- 
tion (blue graph), but it lacks the oscillations evident in 
the quantum curve. They are produced by interference 
among waves traveling along various classical paths from 
the source to each point on the detector. A primitive 
semiclassical calculation (Section llVjl 33l l3^ describes 
the interference (center column). In this calculation the 
wave function at each point is a sum, in which every clas- 
sical trajectory linking the source to the specified point 
contributes one term: Its magnitude is the square root of 
the classical density associated with that trajectory and 
its neighbors, and its phase is given by the classical ac- 
tion (plus a correction — ^ , where p denotes the Maslov 
index for that path). An interesting feature of the semi- 
classical distribution is its approximate symmetry with 
respect to closed orbits (green and purple curves). Hence, 
their significance extends beyond their role in the total 
photodetachment cross section |^ 0| . 

A quantum calculation (Section that makes no ref- 
erence to classical orbits is shown in the right column 
of Fig. n For a point source, the wave function is pro- 
portional to the energy Green function G(f ; E) for an 
electron in the external potential. This function can be 
expanded into a rapidly convergent series over Landau 
states, the eigenstates of the cyclotron motion [T^l 135^. 
This expansion yields both the total detachment cross 
section and the spatial electron distribution. Although 
the mathematical structure of the quantum calculation 
bears little resemblance to the classical analysis, the 
semiclassical results approximate the quantum results re- 
markably well, as a comparison between the center and 
right panels in Fig. ^reveals. The only exception is at the 
caustic surfaces where the classical density diverges. Uni- 
form approximations (not shown) correct the divergent 
behavior near caustics, and result in density profiles that 
are practically indistinguishable from the exact quantum 
result. The very satisfying agreement of quantum and 
semiclassical results for this problem is further illustrated 
in Section IVTl Thus the features of the quantum density 



distribution are best explained from the classical trajec- 
tory fields. 

The purpose of this paper is to present a complete 
analysis of the classical trajectories and their associated 
caustics, and to use that analysis to compute predictions 
of the features that might be seen in a photodetachment 
microscope in parallel electric and magnetic fields. The 
equation of motion is trivial: A superposition of ballis- 
tic motion in the z-direction, and cyclotron motion in 
the X — y plane. Complex behavior only arises when we 
consider the dynamics of the whole family of electron 
trajectories, most naturally formulated as an exercise in 
Hamilton-Jacobi theory. The hardest part of the prob- 
lem is understanding the structure of the caustics — they 
are either delightfully or painfully complex, and a long 
analysis is required to describe them quantitatively. Once 
this is done, it is easy to construct primitive semiclas- 
sical approximations for the wave function everywhere, 
and to derive uniform approximations that correct the 
divergences near the caustics for p > 0. 

In addition, we show how to compute a quantum wave 
function for this system, and we demonstrate that the 
semiclassical approximations are usually in very good 
agreement with the quantum solution. 

II. CLASSICAL MECHANICS 

We first review the properties of the classical motion of 
an electron (charge q — — e, mass m) that is accelerated 
toward positive z in the fields £ = —Se^ and B = —Be^, 
and choose as corresponding electromagnetic potentials 
(<i>,A): 

$(f) = 5z, A(f) = i(Bxf) = |(y,-i,0)^. (1) 

(Physical coordinates f = {x,y,z) are marked with ac- 
cents to distinguish them from their scaled dimension- 
less counterparts introduced below. Primed symbols 
f' = {x',y',z') refer to the source rather than the des- 
tination.) 

The electrons are emitted from the source with fixed 
positive kinetic energy E, i.e., constant initial velocity 
vq = y^2E/m, but arbitrary initial direction, defined by 
their polar angles 9' and (j) ' . ( In the quantum case, tun- 
neling sources with E < are also admissible.) The 
spiraling motion of the electron in the magnetic field pro- 
vides natural frequency and length scales in the problem, 
the Larmor frequency ujl — eB/2m and the maximum 
cyclotron orbit diameter d — vq/ujl, so we introduce di- 
mensionless time and position variables via: 

. t X y z P (o\ 

t = ULt, 2;=-, y=2' ^^d' ^^d' 

In this fashion, the dependence of the result on the pa- 
rameters is eliminated except for one remaining variable, 
which we choose to be the ratio rj of the magnetic and 
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electric forces in the system: 

77 = vaB/E . 

A. Equations of motion 



(3) 



We first characterize the trajectories in the problem. 
The standard minimal coupling Lagrangian £(f , df/dt) 
with the potentials Q is 



m / df \ df ^ 

^=2 U -'Z^W + .^-AW 



2^" 



2 

— z + ■ 



xy 



(4) 



(where the dots indicate differentiation with respect to t) . 
We introduce a complex coordinate s = /5e"^ = x + iy and 
its scaled equivalent s = s/d. The equations of motion 
then read: 



+ 2is = 0, z-2/77 = 0, 



(5) 



and their integration is straightforward. For a particle 
with emission angles (6' , <j)'), we find 



s{t) = sin6''sinte'( 
z{t) ^ f /rj + t cose' + z' 



(6) 



The second of these equations is the familiar uniform 
acceleration in z-direction, while the first describes the 
circular cyclotron orbits in the magnetic field. If the 
source is located at the origin (s' — 0), the motion may 
be expressed in cylindrical coordinates p and (p instead: 

p{t) = smO' \s\i\t\ , (j3{t) ^ (j)' - {t mod tt) . (7) 

The polar angle 4'{t) changes linearly during each cy- 
clotron orbit cycle kir < t < (A: + l)7r, but discontinuously 
jumps by tt at each return to the z-axis (t = kn) ts^]. 
Differentiation of Eq. JHJl yields the particle velocity: 



s{t) = sinei'e'^-^'-^*) , = 2t/ri + cos 6' 



(8) 



Again, the polar angle of the velocity s(t) changes uni- 
formly, albeit with twice the frequency. Finally, we 
also state the (scaled) canonical momentum components 
Px = dC/di = voPx etc. connected with these trajecto- 
ries. From Eq. we find: 



p^[t) = mvQ (s + is) = mt;Q(sin6'' coste"^"^'"*^ -|- is') 
Pz{t) = mvQZ = mvQ {2t/ri + cos 9') . 



(9) 



Action functionals 



The equations of motion (|SJ| describe the path of a 
charge as a function of its emission angles {9' , (f)') and ini- 
tial velocity wq, and the temporal variable t parametrizes 



the trajectory. However, in an experiment, the electron 
travels from a known initial point f' to an observed fi- 
nal point f on the detector. Hamilton-Jacobi theory is 
especially suitable dealing with paths defined by such ini- 
tial and final conditions, and an appropriate quantity to 
carry out calculations is the action integral iS(f , f'; T) for 
a general source point r' and field point r. 

By definition, the action for a given path is the inte- 
gral over the Lagrangian £(f,df/dt) I^J. Inserting the 
trajectory ©, we find initially: 



AiC{Y, df/dt) 



mu^/ST^ AT{T cos 9' + z') 
2a;L V 377^ 77 



Tcos 



2 Ql 



(10) 



sin T cos T sin^ 9' + 2{x'y - xy') 



where T = uj^T. Next, we must eliminate the emission 
angle 9' by demanding that the trajectory r{t) ^ passes 
through the destination f after a time of flight T — for 
instance, from Eq. Q we find sin^ 9' = {s - s'Y / sit? T. 
The action then can be written as the sum of a longitu- 
dinal and a transverse part, 



5(r,r';T)=5i(s,s';T)+5ii(z,z';T) 



(11) 



which read respectively: 



5j_(s, s'; T)='p^ (|s - sf cot T + 2{x'y - xy')) 
2uJL 



(12) 



S\\{z,z'-T) 



2lul 



{z~z')^ 
T 



2{z + z')T 

rj irf 



(13) 

These results for the classical action in a two-dimensional 
magnetic field, and for the linear potential in one dimen- 
sion, are well known |34]. 

Because 5(r, r';T) is a single- valued function, exactly 
one trajectory with given time of flight T will link r and 
r'. Its initial momentum p' = —V S is given by: 



f OS .dS\ 2r/ /^ m • 1 



\dx 
dS_ 



T 



(14) 



We now assume r' = 0. The initial momentum along 
the z-axis is p'^ = mv^ cos9' , so the emission angle as a 
function of r and T follows from Eqs. Q)-©: 



^'(r; T) = arccos(z/T - T/rj) , 
(j)'{r;T) = (f>+ {T mod tt) , 



(15) 



where </> is the polar angle of the destination r. Similarly, 
the flnal momentum p of the particle at its destination r 
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is given by the derivative V5: 

fdS .dS\ i\ m • /n 

" '^^ V ' 9^ J " """"o [(s - s ) cot T + IS J , 



Pz = t^L 



95 

dz 



T 



(16) 



In particular, for r' = the velocity component in z- 
direction at r @ reads: 



i(r;T) = z/T + T/r^ 



(17) 



We define E{y, f '; T) as the energy of the electron that 
goes from f' to f in time T. Similarly, E{y\ T) is the 
energy of a particle that starts from the origin (f ' = 0) 
and arrives at the scaled position r at the scaled time 
T . These energy expressions follow from the derivative 
of the classical action: E{r,r';f) = -dS{r,r';f)/df, so 
Eqs. I)ll|l -H13 |l yield the decomposition: 



with 



E«{z-T) 



sin^T 



(18) 

(19) 
(20) 



The experiment provides a stationary source that emits 
particles with fixed energy E = mwg/2, rather than with 
a definite time of flight T. Therefore, only those classi- 
cal trajectories fulfilling the condition E = E{r; T) are 
acceptable solutions. Hence, we may interpret Eq. lfTH|) 
as an implicit equation for the time of flight T{y\ E) as a 
function of the particle energy: 



6(r) 



E{v,T) 
E 



sin^T Vr 



T\2 



(21) 



While this transcendental equation must be solved nu- 
merically for T, this is the only instance in the classical 
description of motion in parallel fields that we have to 
resort to computation. Every other quantity of relevance 
can be cast as a closed-form expression involving the time 
of flight T. For example, the reduced action (Hamilton's 
characteristic function) W{r; E) that governs the semi- 
classical phase follows from Eq. 1)21(1 together with the 
time-dependent action 111(1 via: 



W{v-E)^ f p-dq = S{v;f{r;E))+Efir;E) 
Jo 



(22) 



In general, Eq. possesses a variety of solutions T. 
Correspondingly, the reduced action yV{r;E), like other 
quantities of interest in the stationary theory, is a multi- 
valued function and depends on the specific trajectory 
chosen. 



C. Number of trajectories 

Clearly, the multiplicity of real solutions T to the im- 
plicit equation 1(21(1 determines the number of classical 
trajectories with given energy E that connect r' = 
with r in parallel fields. The semiclassical analysis of 
the problem rests on our knowledge of these trajectories. 
Here, we establish some general properties of the set of 
solutions, js^l 

To learn more about Eq. ((21(1 is is useful to sketch the 
scaled particle energy e{T) as a function of the time of 
fiight T (Fig. [2J|. In the graph, we plot the longitudi- 
nal contribution to the energy e||(r) l(2U(l (red line). It 
approaches infinity both for T —^ and T — > oo, runs 
for z > smoothly through a minimum e(Tff) = at 
Tff = ^/rfz which corresponds to the "freely falling" mo- 
tion of a charge starting at rest, and intersects the line 
£||(T±) = 1 at the scaled times T±: 



T+ 



yrf^+Arjz ± rj 



(23) 



These are the flight times of electrons emitted parallel to 
the flelds. For z>Q their difference AT = 77 is constant 
and equal to the return time for the parallel closed orbit. 
(Note that in physical coordinates, T± ~ T±/lol docs 
not depend on the magnetic fleld.) Since e(T) > e||(T), 
all real solutions of Eq. ((21(1 are located in the interval 
T_ < T < T+. 

Adding the perpendicular component ej_(r), Eq. 1(19(1 . 
yields the full dependence of e(T) on T (black curves 
in Fig. 121). For p > 0, this "magnetic" part of the 
energy diverges whenever the time of flight Tk = kir 
[k = 1, 2, 3, . . .), is a multiple of the cyclotron orbit time, 
when every trajectory must return to the z-axis. Halfway 
between these singularities, ei(T) adopts its minimum 
value e±_ [(fc — ^)7r] — ■ For p> 1, no real solutions to 
Eq. ((21(1 exist, and the sector of classically allowed motion 
is entirely contained inside (but not filling) the cylinder 
p < 1. It is a simple exercise to verify that within each 
cyclotron orbit cycle [k — \)-k < T < kn, e(T) is a con- 
vex function that has a single minimum e^^j^j^. Its value 
determines the number of trajectories in this cycle: For 



mill 



< 1 there exist exactly two times of flight T^*^-* and 
Ty'"' that correspond to an early and a late trajectory 
that both complete k — 1 cyclotron orbits before arriving 
at r (marked by x in Fig. |2J|. 

With increasing radial distance p, e'^J^^^^ grows until it 
becomes equal to unity, causing the two trajectories to 
coalesce into a single solution. For even greater distances 
p, e|^;^jj > 1 holds, and no classical trajectories exist in 
the cycle. (This transition spawns a complex conjugate 
pair of solutions that make up a tunneling trajectory.) 
We infer that with increasing radius p, the number of 
solutions N{p, z) successively drops by two. (See Fig.^) 

Consequently, for any given distance z, the maximum 
number of trajectories occurs in the vicinity of the z-axis 
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FIG. 2: Graphical solution for the number of trajectories from 
source to detector. The heavy (red) curve is t\\{T), Eq. (1191 . 
Its intersections with the line e = 1 define T- and T+ (circles), 
the times of flight of electrons respectively emitted in the 
direction, and against the direction, of the electric force. The 
light (black) curves trace e{T), Eq. (HTJ. Each goes to infinity 
at T = kn, i.e., whenever a cyclotron orbit is complete. Every 
intersection of these curves with e — 1 (crosses) represents the 
time of flight of a trajectory from the source to the specified 
point r on the detector. There are at most two solutions 
within each cyclotron period (fc — l)7r<T<fc7r. (In this 
graph, 77 = 11, 2 = 3, and p — 0.3.) 



(p — > 0). For infinitesimal p, the lateral energy e±{T) be- 
comes negligible, except in the immediate neighborhood 
of the return times Tk, where e±{T) still diverges. As a 
result, the late and early trajectories Ty~^\ T^^ in con- 
secutive cycles generally tend towards the singularities 
at Tk = kn, while the fastest and the slowest solutions 
approach the parallel times T± (|23|l . respectively: For 
/9 — > 0, every cyclotron orbit period (fc — l)7r < T < kir 
that intersects the interval < T < contributes 
two real solutions to Eq. 121|) . Thus, the total number of 
solutions N{p 0,2) depends only on the ratio of the 
interval sizes, and their relative positions. It is easy to 
see that for z > 0, 



iV(p^0,z) = 2[Lr/AJ-fn(z)] 



(24) 



where [. . .J denotes the integer part, and n{z) is cither 
1 or 2, depending on the distance z. (In particular, for 
Tft = kn, n{z) = 2 holds, while niz) = 1 for Ts = 
{k — ^)n.) We infer that the ratio of forces ry presents a 
convenient control parameter that adjusts the degeneracy 
of the solution set. We also note that no matter how small 
we choose the magnetic field B, in the vicinity of the 
"focal points" z — fc^Tr^/ry always a region with at least 
four coexisting trajectories exists (23] ■ FigureHshows an 
example with up to eight interfering trajectories. 

For a destination on the z-axis (p = 0), the solutions 



In this case, the polar angle of the destination (p is no 
longer defined, and the emission angle (j)' (|15|l is unre- 
stricted: An entire cone of trajectories connects the origin 
with r, reflecting the cylindrical symmetry of the situa- 
tion. In particular, this covers the situation r = 0, where 
the theory yields the orbits returning to the source: Be- 
sides the "uphill" closed orbit emitted antiparallel to the 
electric force, having time of flight AT — rj H23|) and 
action >Voi(0;S) = 2r]E/{3ujL) (EU, there exist [rj/nl 
additional "magnetic" closed orbits with time of flight 
Tfe — kn, emission angle 9'f. — arccos(— fcTr/yy), and ac- 
tion: 



>viln(0;i?) 



knE 



(25) 



T 



(fe-i) 



and T^F' coincide with the orbit time Tk = kn. 



where k < ^jn. In (p, z) space, these orbits have the 
shape of either "balloons" (even k) or "snakes" (odd fc). 
For the latter, both p Q and z (O are zero at i = T'fc/2: 
The electron stops, turns around, and retraces its path 
back to the source. (See also Fig. Q]) A semiclassical 
study of the total photodetachment cross section in par- 
allel fields has been based on these returning trajectories 



III. THE CAUSTIC SURFACE 

Placing the detector at fixed distance z from the 
source, as we move outward in p, at certain discrete dis- 
tances the number of classical trajectories arriving at the 
point (p, z) drops. These locations, at which an infinites- 
imal shift changes the number of solutions, make up the 
boundaries of classically allowed motion, the caustic sur- 
faces js^. A simple way of illustrating the caustics is 
depicted in Fig. ^ (top left panel): Using the equation of 
motion (O, |7J, we draw a series of classical trajectories 
with increasing emission angle d' (blue curves). These 
trajectories then trace out their bounding surfaces (red 
lines). The result is a structure of considerable complex- 
ity. (Similar presentations, for varying values of the ratio 
of forces 77 can be found in Ref. ^.) 

We would like to obtain a closed functional form p(z) 
representing the caustics, but no such explicit depen- 
dence can be found. Instead, we employ a temporal 
parametrization of the caustic surfaces [p(t),z(t)] which 
yields the location of a point on the caustic as a func- 
tion of the time of flight t required to reach it from the 
source. In this parametrization, each cyclotron period 
(k — \)n < t < kn gives rise to an individual caustic 
surface. Additional focal line segments on the symme- 
try axis (where trajectories cross through p = after 
k complete cyclotron periods, i.e., t — kn) also may be 
regarded as caustics. 

The resulting caustic surfaces themselves have a sur- 
prising degree of complexity: Seven different types oc- 
cur (Fig. ^ that regularly change into each other with 
increasing parameter rj. This evolution of the caustic 
structure is illustrated in Fig. O These results follow 
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from an elementary but lengthy argument based on the 
irregular points of the parametrization (that may or may 
not show up as singularities in the caustic curve in space 
p{z)). We only sketch it here; a complete proof is given 
in the appendix. The reader who is not interested in 
mathematical detail may simply scrutinize Figs. 0] and El 
together with Tabled and then skip to Section Hvl 



A. Temporal parametrization of caustics 



where: 



tani 



9 

77^ r 



A{t, t)^T^ --'—+7] 



(30) 



Signs must match in the pair of equations H29|l . The 
corresponding emission angle follows from Eq. IjlSf) : 



cos 6'^{t) 



i(r± V^(t, t) 
ri{t - t) 



(31) 



A point on the caustic surface r is characterized by two 
coalescent trajectories connecting it to the source, which 
we choose to be the origin. In the time-of-flight picture, 
this requires not only that the energy £'(r; t) H18|l for 
the trajectory matches the particle energy (e(t) — 1), 
but that additionally E{t; t) as a function of t attains a 
minimum there (two roots of e{t) = 1 coalesce). Hence, 
f will be located on the caustic surface if: 



at 



and 



^(f;f)^o 



hold simultaneously, or equivalently, 

de 

e{t) = 1 and ^ = ■ 



(26) 



(27) 



We already examined the first requirement, Eq. (|21|) . in 
Section Hi CI From Eqs. lfTT|) - ({T^ . the second condition 
reads in explicit form: 



1 de cos t 
'2dt~ sin^t 



z2 t 



t3 



(28) 



The set of Eqs. and for the three variables 
p, z, and the time of flight t implicitly characterize the 
caustic surfaces in the parallel field problem. It would 
be nice if we could eliminate t between these equations, 
and represent the caustics more directly as Pcaustic(-2)- 
(Indeed, this strategy works for a uniform pure electric 
field; one finds that the caustic is a paraboloid with the 
source as its focus |l7j.'l However, for nonzero magnetic 
field, the trigonometric quantities in Eqs. (|21() and H28() 
prevent elimination of t. 

Instead, we can get a parametric representation of the 
caustics by the following method. Using Eq. we 
express as a function of z and t, and then insert the 
result into Eq. H21|l : this yields a quadratic equation re- 
lating z and t. We put the solution to this equation back 
into Eq. (|28|l . and obtain a temporal parametrization of 
the caustic surfaces: 



z±{t) 



rj{t - r) 



P±it) = 



I sinfl 



1 - 



t±^A{T,t) 



t - T 



(29) 



B. Existence of solutions 

The expression (|29|l shows that for a given time of 
flight i, there exist at most two corresponding points 
[p±{t), z±{t)^ on the caustic surface. These solutions are 
not necessarily real, however. In order to qualify as a 
physical solution, a pair of inequalities must be fulfilled. 



A{t, t)>0 and {t 



forT>0, 
(32a) 



or 



A{t, t)>0 and {t - r f < {t ± ^/^(t,*)) /or t < . 

(32b) 

The signs "±" refer to the branches in Eq. H29|) . Note 
that the ordering in the second inequality depends on 
the sign of r = tani. (For r = 0, both solutions exist 
and lie on the symmetry axis p = 0.) 

An algorithm for selecting the upper (+) and/or lower 
(— ) sign in Eqs. consistent with these condi- 

tions is given in Fig. |21 To understand it, it is best to 
think of A{t, t) H3Q|I as a function of two independent 
variables. For t > f' ^('''i^) never negative, but for 
t < ^, A(T,t) usually has two roots in each cyclotron 
period (k — 1)tt < t < kn, the larger of which is: 



4i2 



(33) 



(The period containing i = is the only possible ex- 
ception.) We call the solutions of this transcendental 
equation the minimal times: 



Jfc) 



(34) 



With this information, the algorithm is easy to apply. A 
proof is given in the appendix. 

As t increases from to |, r = tani sweeps from 
toward infinity, and then from — cxd to +00 an integer 
number of times, finally ending up (usually) at some finite 
value of T. In this range of t, whenever r lies between 
and Ty \ both (+) and (— ) solutions are complex. Two 

(k) 

real solutions are created at , and they persist while 
T increases to -l-oo, jumps back to —00 and increases 
again, until r passes through zero, at which point they 
disappear (Fig. EJ. 
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Compare t, r|/2 



f<Tl/2 




t=r[/2 




f>r|/2 



< tan(ri/2) < rj 













no 




yes 


















tan f < 




tan t= 




tan f > 



T 






3 < tan f < x'^>(0 ?\^ 













no 




yes 







X 



FIG. 3: An algorithm detailing the existence of real solutions 
in Eq. 



Once t exceeds ^, either the (+) or the (— ) solution is 
real, according to whether r < or r > 0, respectively, 
so p{t\ z[t) are single-valued functions of i. (+) and (— ) 
solutions join each other smoothly when r passes through 
infinity. 



C. Irregular points on caustics 

Usually, the parametrization [p±(i), z±(t)] H29I) is a 
smooth function of the time of fiight <, and the resulting 
caustic surface in space is locally smooth as well. We call 
such points on the caustic regular. For some isolated val- 
ues to, however, one or both of the functions p±{t), z±{t) 
become singular. While such irregular points (io,2^0jPo) 
are rare, they are important for our understanding of 



the caustic structure because qualitative changes of the 
parametrization may only occur at these points. Like- 
wise, singularities of the caustic surfaces themselves are 
always connected to irregular points. 

Because they are involved with changes in the num- 
ber or the branch of solutions, the irregular points of 
the parametrization can be read off the flowchart algo- 
rithm in Fig. O We identify eight different types: (Case 
designations refer to the corresponding section in the ap- 
pendix; the symbols are used in the figures.) 

a. Maxima of p{t) (□). The irregular points are at 
to ^ {k~ ^)Tr, where r — tani diverges and changes sign. 
Here, the (+) and (— ) solutions switch seamlessly into 
each other. On the caustic surface, the corresponding 
irregular points mark the maxima of p (case 2a) , as well 
as the endpoints of snake- type closed orbits (case 2b). 
The latter exist only for < < ^. 

b. Top of caustic (★). This irregular point is at 
(toj^OiPo) — (f'^f'*^)' uphill turning point (case 
2c). There, the radicand in p±(t) vanishes for either the 
(-I-) or the (— ) solution (the remaining solution is unaf- 
fected), and the number of physical solutions changes by 
one. The caustic surface itself is usually smooth. 

c. Minimal times (0). Here, to = t^^^^ ^'M^ is a root 
of A{tant,t) = (case 2d). The singularity occurs only 
for t < ^ and indicates a minimum of the time of fiight 
from the source to the caustic. At t = to, the solution 
pair [p±{t) , z±{t)] springs into existence. 

d. Complete cyclotron orbits (o). For to = kir the 
caustics intersect the symmetry axis {po — 0), and tant 
changes sign (case 2e). Caustic surfaces "begin" and 
"end" here. In contrast to the previous cases, the irreg- 
ular points at zq = ^{kn ± r]) mark singularities in the 
caustic profile p{z). Usually, the caustic surface locally 
has the shape of a rotationally symmetric cylindrical cusp 

Focal line segments, parts of the caustic structure 
where the trajectories are focused back onto the symme- 
try axis, connect the cusp singularities. 

These irregular points are "generic," i.e., they occur for 
almost all choices of the parameter 77. At some isolated 
values of 77, however, the top of the caustic {to = §) 
coincides with another irregular point, and the ensuing 
singularity has distinct properties. Such values of rj are 
important when we try to figure out how the pattern of 
caustics evolves with increasing r]. Three cases arise: 

e. Hyperbolic umbilic points (x). When rj = 2fc7r 
and to = kn (case 3a), the time to reach the top of 
the caustic is a multiple of the cyclotron period, and the 
uphill turning point is shared by two locally degenerate 
caustic surfaces in the shape of a cylindrical cone with 
opening angle j. They are manifestations of a higher- 
order singularity known as hyperbolic umbilic points |37| . 
Also, a "balloon" -type closed orbit is created as 77 in- 
creases through 2k7T 9]. 

/. Bifurcation points (v)- Here, to — t^^^^ — ^ holds 
for the minimal time (case 3b). Then, the force parame- 
ter rj takes on one of its critical values rj^^^^^ given by the 
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roots of the transcendental equation: 

Vi,, = tan{r^ij2) , (35) 

where k — 1,2,... refers to the interval 2(fc — l)7r < 
7y^j.it < {2k — l)n. At these points, a bifurcation occurs 
that affects the temporal ordering of the uppermost caus- 
tic, but not its shape. 

g. Creation of a closed orbit A. Finally, for rj = 
{2k — l)7r (case 3c), the time of flight for the antiparallel 
"uphill" closed orbit is an odd multiple of the cyclotron 
period. At this value of rj, a new "snake" -type closed 
orbit is created 9J. 

To obtain these properties, we studied the local behav- 
ior of the parametrization (|29|l near the irregular points. 
Details are discussed in the appendix. 



D. Classification of caustics 

We listed eight different types of irregular points in the 
parametrization of the caustics. However, when we look 
at the caustic surfaces in (p, z)-space, we find, happily, 
only two shapes of caustics that are present for general 
values of rj, and a third occurring when 77 is a multiple of 
2tt. They all resemble Christmas tree ornaments cusped 
at their bottom, but differing at their tops: The upper- 
most region may be a cylindrical cusp shaped like an 
onion-dome (type [O]; the slope dz{p)/dp goes to infin- 
ity), or a conical lid (type [C]; the slope approaches unit 
value at the uppermost point), or a smooth dome (or in- 
verted tear-drop, type [S], with zero slope). The first two 
basic shapes come in two variants each: In the "early" 
(e) or "closed-orbit" version, a focal line runs through the 
centerline of the ornament {p = 0) all the way from the 
bottom to the top. They also harbor the endpoint of a 
"snake" closed orbit. In the "late" (1) or no-closed-orbit 
variant, the focal line goes only part way through the 
ornament, from the cusp at its bottom to a point within 
the ornament where it connects to an upward-pointing 
cusp from another ornament. The three teardrop-shaped 
caustic subtypes share the latter property; they differ 
only in their parametrization. All seven types are shown 
in Fig. 21 



1. Onion-domes 

a. Late variant [01]. The caustic created in the fcth 
cyclotron orbit ((fc — l)7r < t < kn) belongs to this class 
if i > ^ holds for all points of the caustic. In other 
words, the fcth caustic surface is a late onion-dome if 
77 < 2(fc — l)7r. Its parametrization is simple: The caustic 
starts at <o = (^ ^ l)7r from its upper cusp at po — 0, 
zo = ^0(^0 ~ il)/v (case 2e). As t increases, p-{t) and 
z_(t) trace out the upper surface of the "onion." When 
t approaches {k — 5)77, P-{t) tends towards its maximum 
value p = 1 (case 2a). The (— ) solution ceases to exist 



there, and the (-I-) solution smoothly takes over. It then 
traces out the lower part of the surface down to the lower 
cusp point to = kn, po — Q, zq — ^{kir + r/) (case 2e). 
From there, a focal line of length 2fc7r runs back along 
the symmetry axis to the upper cusp zq = ^{kn — if) 
of the next late onion-dome, creating an infinite chain of 
overlapping surfaces linked by focal lines. 

b. Early variant [Oe]. If t < holds for all points of 
the fcth caustic surface, i.e., 77 > 2fc7r, the caustic has the 
same shape, but a more complex temporal parametriza- 
tion: We begin at = kn with the (— ) solution at the 
top cusp of the onion at zq ~ ^{kix — r\). However, we 

now allow t to decrease. When t reaches (fc — ■^)tt, we 
have arrived at the endpoint of a "snake" closed orbit 
that returns to the origin (case 2b). To follow the curve 
smoothly, we switch to the (+) solution, and then con- 

(k) 

tinue downward to to — ^,,-^4, the minimal time of flight 
to the caustic (case 2d). There we switch back to the (— ) 
solution, and now allow t to increase. When t approaches 
{k — i)7r, we are at the maximum radial extension p = 1 
(case 2a). We switch again to the (-I-) solution which 
describes the entire lower portion of the caustic down to 
the bottom cusp at zq = ^{kn + rf) that is attained 
for to = kn (case 2e). Note that the parametrization 
[/9(t),z(t)] for this type is a double-valued function of t. 

Again, a focal line of length 2fc7r runs along the center 
of the ornament, but in this case it connects the opposite 
cusps of the caustic itself. Early caustics, in contrast 
to their late counterparts, are isolated structures — no 
continuous curve {t,z,p) on the caustic set will connect 
them to another caustic surface. 



2. Conical lids 

Apart from the cone-like shape of the caustic at the 
upper end point (case 3a), they share all other aspects 
with the onion-dome caustics discussed above. Conical 
lids occur only when rj is an integer multiple of 27r. For 
rj = 2fc7r, the (fc -|- l)th caustic surface turns into the late 
variant of the conical lid [CI] that has the pattern of a late 
onion-dome, and is linked to their sequence via a focal 
line segment. Simultaneously, the fcth caustic becomes 
a conical lid in its early variant [Ce], a self-contained 
caustic with an associated closed orbit that is patterned 
after the early onion-domes. The caustics are joined at 
the tip of the umbilic double cone structure. A close-up 
is shown in Fig. 1111 (Appendix) . 

3. Smooth domes (inverted tear-drops) 

For 2(fc — 1)tt < rj < 2fc7r, the fcth turning surface 
(with times of flight between (fc — l)7r and kir) becomes a 
smooth-dome caustic that forms the uppermost bound- 
ary of classical motion, terminating on the symmetry axis 
for to = 2 Zq = — J. We may discern three types that 
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(k-l)jt tan/ k7T (k-1)7i: tan f kTU (k-1}7c tan f kTi (k-1}7i: tan f kjc 



[01] : Onion-Dome [Oe] : Onion-Dome [CI] : Conical Lid [Ce] : Conical Lid 

(late variant) (early variant) (late variant) (early variant) 




[SI]: Smooth Dome [Si]: Smooth Dome [Se]: Smooth Dome 

(late variant) (intermediate variant) (early variant) 



FIG. 4: The seven types of caustic surfaces occurring in the parallel field problem. They form three families with distinctive 
shapes of the turning surface: Omon-domes (top left pair) come in late (or no-closed-orbit) variants (left, rj < 2(k — and 
early (or closed-orbit) variants (right, r) > 2kn). The related late and early conical ltds (top right pair) are encountered for 
rj = 2(k — l)n and rj = 2k-K, respectively. The bottom row depicts the smooth domes (or inverted teardrops) occurring for 
2(k — < rj < 2kTT. The three subtypes differ only in their parametrization, as explained in the text, (fe denotes the index 
of the caustic; i.e., electrons arrive at the caustic in their fcth cyclotron orbit, (fc — l)n < t < kir.) — Beneath the picture of 
each caustic is a time axis showing the temporal parametrization l|29|l . As t follows the green curve and we switch between the 
(-I-) solution (red) and the (— ) solution (blue), the path in space traces the caustic from the highest point to the lower cusp 
(green arrows). Orange parts of the symmetry axis represent focal lines of infinite degeneracy [t — kn), and purple curves are 
closed orbits. Symbols represent irregular points of the parametrization (see Appendix): Half-integer orbits (□, cases 2a, 2b), 
top of caustic (*, case 2c), the minimal times t^^^^ (0, case 2d), cylindrical cusps (o, case 2e), and umbilical points (x, case 
3a). Black dots (•) mark the source location. — Parameters used (left to right): rj = 1.5, k = 2; rj = 7.908, k — 1 (cf. Fig. 0; 
r] — 2-K, k = 2 and A; = 1 (cf. Fig. Illll (top row); r/ = 1.5, 77 = 2.75, 77 = 4, all fc = 1 (bottom row). (Force directed downward, 
caustics not to scale.) 
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interpolate between the early and late regimes of the mo- 
tion (Fig.^ bottom row). In each case, a focal line seg- 
ment joins the lower cusp of the caustic to the uppermost 
cusp of a late onion-dome surface. 

a. Late variant [SI]. This case holds for 2{k— l)n < 
V < Vait ISSll- We start with the (— ) solution from the 
top of the dome (fg = f ' '^^^^ 2*^)- ^ increases, we trace 
out the upper surface, and when t goes through (A; — i)7r, 
and p reaches its maximum (case 2a), we smoothly switch 
to the (-I-) solution that traces out the lower surface. The 
caustic curve is a single- valued function of t. 

b. Intermediate variant [Si]. Here, ry^^j^ < rj < [2k~ 
l)7r. We now start at io = f fro™^ the uppermost point 
with the (-I-) solution. As t decreases to t^^^.^ (case 2d) we 
trace out part of the smooth dome. Then we switch to the 
(— ) solution and allow t to increase again to (fc— i)7r (case 
2a), completing the upper half of the caustic. Switching 
back to the (-f ) solution we trace the lower part of the 
ornament down to the cusp for t^ = kir. 

c. Early variant [Se]. In this case we again begin 
with the (— ) solution at the dome top {to = |), but now 
{2k — l)7r < rj < 2kTr holds. As t decreases through (fc — 
^)7r (case 2b), we switch smoothly from the (— ) branch to 
the (-I-) branch at the endpoint of a snake-shaped closed 
orbit. The time t continues to decrease, and construction 
of the remaining part of the caustic follows the pattern 
of the intermediate case sketched above. 



E. Evolution of the caustic pattern 

Considering the simplicity of the trajectories, the caus- 
tic surfaces in the parallel field problem are of surprising 
complexity. Let us now examine how the caustic pattern 
evolves as we increase the energy of the emitted electron, 
or, more generally, increase the parameter rj = vqB/£. 
Starting with « 0, we obtain the sequence of changes 
displayed in Fig. 

a. < 77 < rjl^i^. {Vait ~ 2.331 is defined in 
Eq. 1)35(1 .) The top of the first caustic {k = 1) is a smooth 
dome, and the curve is parametrized in the way of the 
late variant [SI]. Its lowest point is a downhill-pointing 
cusp corresponding to io = at zq = ^{n + r]) (Eq. lO). 
Extending upward from this cusp is a focal line, which 
ends at Zq = ^{ti — rj). This is the uppermost point 
of an onion-dome caustic that begins here in an uphill- 
pointing cusp. Its parametrization is again the late vari- 
ant [01]. Below it hang an infinite sequence of similar 
onion-domes interconnected by focal lines running from 
the lowest point of each caustic at zq = ^{kT^ + v) to the 
highest point of the next one at zq = ^{kn — r/). Near 
these focal lines, four trajectories go through every point 
on the detector. 

h. rj = We are at the first bifurcation 

point (case 3b). The shape stays the same, but the 
parametrization of the smooth dome changes here. 
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TABLE L Summary of properties of the caustics listed in 
Fig. 0] The "early" or "closed-orbit" types [Se], [Ce], and 
[Oe] contain the endpoint of a snake-type closed orbit. Early 
onion-domes [Oe] are additionally associated with a balloon- 
type closed orbit. This type is also isolated; the focal line 
associated with it joins its two cusp singularities on the sym- 
metry axis instead of connecting the lower cusp to the oppo- 
site cusp of the following "late" onion-dome [Ol] . (The conical 
lids [CI], [Ce] marked by * are linked at the umbilical point.) 
The late types [01], [CI], [SI], together with the intermediate 
case [Si], all share downhill motion, i.e., z{t) > 1171 1 holds 
for all points of the caustic. The uphill turning point belongs 
either to one of the smooth dome caustics [Sl]-[Se] or to the 
pair of conical lid types [CI], [Ce]. For the types [Si] through 
[Oe] , the minimum time of flight to the caustic occurs off the 
symmetry axis, while the late types [01]-[S1] are single- valued 
functions of t. — The last row illustrates that the caustic 
types in the given order form an evolutionary sequence: With 
increasing 77, each caustic is transformed in six successive steps 
from a late onion-dome [Ol] to an early onion-dome [Oe]. 

c. rjl^-it < 7] < TT. The parametrization of the 
smooth dome is now the intermediate-time variant [Si], 
with a newly created off-axis minimum (case 2d) in the 
time of flight. As rj increases, the lengths of the onion- 
dome caustics increase, with their upper cusps moving 
up and their lower cusps moving down. 

d. rj = n. Several things happen here. The top of 
the highest onion-dome moves up to the origin, and a 
snake- type closed orbit is created (case 3c). The bottom 
cusp of each ornament touches the cusp at the top of the 
second ornament down. Therefore, focal lines now cover 
the entire positive z axis. The parametrization of the 
smooth dome changes again. 

e. TT < rj < 2tt. The uppermost caustic is now 
parametrized by the early variant [Se] of the smooth 
dome. The first snake-type closed orbit is present, and 
the cusp that had touched the source has moved further 
uphill. New diamond-shaped regions appear in which six 
trajectories lead to every point. In these regions, the fo- 
cal line at p = may be said to be degenerate: There 
are two distinct families of trajectories that focus at the 
symmetry axis. 

/. rj = 2tt. We have reached the umbilical point 
(case 3a). The uppermost cusp now has caught up with 
the uppermost point of the smooth dome. The caustics 
involved in the "collision" have been deformed into coni- 
cal lids — the former early-type smooth dome has turned 
into its early variant [Ce] , while the late-type onion dome 
takes on the late- type parametrization [CI] . Another cusp 
has moved up to touch the origin, and bottom and top 
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Radial Distance ^ [|im] 



FIG. 5: Evolution of the caustic surface in parallel fields £ — 15V/m and B — 0.02 T. Left to right: Emission energy 
E = 3^eV, 8.690 /ieV, 12 ^eV, 15.78 peV, 40/ieV, and 63.13 /xeV, corresponding to r? = 1.370, rj^it. ~ 2.331, 2.739, tt, 5.001, 
and 2n, respectively. — Red and blue segments denote the (+) and (— ) solution in the parametrization 12911 . while focal lines 
connecting the caustics are shown in orange. Purple (dashed) curves: Snake-type closed orbits. Symbols indicate irregular 
points of the parametrization: Half-integer orbits (□, cases 2a, 2b), top of caustic (★, case 2c), minimal time of flight (0, case 
2d), cylindrical cusps (o, case 2e), umbilical points (x, case 3a), bifurcation point (v, case 3b), and creation of a closed orbit 
(A, case 3c). The dots (•) mark the source location. Yellow (blue, purple) shaded areas represent destination points connected 
to the source by TV = 2 (4, 6) classical trajectories. Caustics are labeled according to their classiflcation (green). — (The 
caustic surface with E = 100 /leV (77 — 7.908) depicted in Fig. 0would present the next step (g) in the sequence shown.) 



cusps of different ornaments again coincide. This is the 
bifurcation point for a balloon- type closed orbit. 

g. 2tt < T] < ?7^i.it. This situation is depicted in 
Fig. n The former onion dome (fc = 2) has moved to 
the top of the caustic structure and has further evolved 
into a smooth-dome with late-type parametrization [SI], 
while the original smooth dome (k = 1) has become the 
first onion-dome belonging to the early variant [Oe] . The 
ballon-type closed orbit has detached from the symmetry 
axis, and diamond-shaped regions with eightfold path 
degeneracy occur. Apart from the presence of the 
innermost [Oe] caustic, the caustic pattern is the same 
as in step (a). 



ues to evolve in the same manner. Each time 77 increases 
through TT, a new closed orbit is generated, alternating 
snakes (77 = {2k — l)7r) and balloons (77 = 2A:7r), and the 
maximum number of classical paths leading from the ori- 
gin to any off-axis destination (p 7^ 0) increases by two. 
The fcth caustic surface undergoes its evolution in the 
interval 2{k — l)7r < rj < 2k-K that ultimately results in 
the transformation of an onion dome [01] with the late 
parametrization scheme into an onion dome of the early 
variant [Oe] . — The properties of the different types and 
variants of caustics in parallel fields, and their evolution 
with increasing 77, are briefly summarized in Tabled 



At this point, the evolutionary cycle (with a period 27r 
in 77) starts over: As 77 grows further, the pattern contin- 
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IV. SEMICLASSICAL AND UNIFORM 
APPROXIMATIONS 



p = pd ^ and a polar angle cj) |(7J . Particle conservation 
then requires that the emission rate at the source: 



Except at the caustic surfaces, the quantum wave func- 
tion '0(f) generated by a fixed energy point source in par- 
allel fields is well approximated by a sum over all N{r) 
classical paths connecting the source with the destination 
point f on the detector |32|] : 



V'sc(f) 



N 

E 

Q = l 



0a (f) 



(36) 



The magnitude and phase of the semiclassical wave func- 
tion for trajectory a depend on its classical density 

Pa — Pci \^) s-iid classical action functional Wai^; E) 
respectively [s^ : 



V'a(f) = -y/pjFj exp(iXa(f)) 

Here, the phase Xai"^) is given by: 

1 TT 

Xa(f) = -Wair;E) - -Ma(f) 



(37) 



(38) 



The Maslov index Pa{^) keep s track of the singularity 
"history" of the trajectory [33 • A more accurate result is 
obtained if we extend the sum H3f)|l to include tunneling 
trajectories. Their times of flight T(r; E) are complex 
roots of Eq. igll); their action Wa{r]E) formally 
obtained by analytic continuation, has an imaginary part, 
so their contribution to the semiclassical wave function 
is exponentially suppressed. 

The wave function (|36|l leads to semiclassical expres- 
sions for the particle density Psc(f) and current density 
distribution jsc(f). Ignoring tunneling trajectories, they 
differ from their classical counterparts by the presence of 
interference terms |59| : 

N 

Psc(f) ^^Pa+^Y^ ^JPaPfS C0s(Xa - Xp) , (39) 

a=l a<P 
N 

jsc(f) = E + E VPo'PP ('^a + V/3)cOs(Xa - Xp) ■ 

a=l a<f3 

(40) 

Here, Vq, is the electron velocity at the detector f , and 
ia = Pa'^a dcuotes the local classical current density 
along trajectory a. 



A. Classical density and current 

To find the WKB wave function, we now determine the 
particle and current densities pa{r), ja{f) for the classical 
trajectories a in parallel fields. A bundle of trajectories 
emitted from the source under spherical angles 9' and 
4>' H15|) will cross through the detector surface (here a 
perpendicular section with fixed z) at a radial distance 



a. 1 



(41) 



equals the flux through the detector surface (with normal 
vector a): 



J„ = j„(f) • da = j,"(f )pdpd0 = ^j-(f)d(p2)d0 . 

(42) 

Jo denotes the total emission rate which we set equal 
to the total current generated by an isotropic quantum 
point source C(5^(r) of free electrons with wave number 
k = {2mEY''^/h [13: 



ME) = \C\ 



I mk 



(43) 



(See Section rV CI ) f{9',(j)') is the (normalized) angular 
distribution of the source. We assume isotropic emission, 
i.e., an s-wave source with f{9',(j>') — l/(47r). 

Equating (|1T|| . we find for the current density j": 



9l )2 

j-(f) = ^jo/(c,0:.) 

«0 



9(cos0L>L) 



9(p2,0) 



(44) 



The relationship between the emission angle and the loca- 
tion on the detector plane is given by Eq. (|15|) , where we 
have to keep in mind that the time of fiight T = T(r; E) 
in this formula is a function of the distances p and z, and 
the energy E. Because d(j)'/d(j> — 1 and d cos 9' /dcj) — 
the Jacobian in Eq. (|44|l reduces to a simple derivative: 



d cos 9' 



9p2 



(45) 



To evaluate it, we note that (from Eqs. (|12|) and l\lb\i ) 
d cos 6' 



z 
T2 



1\ dT 



z dT 



(46) 



Implicit differentiation of e{T) H21|l with respect to p^ 
yields: 



1 



sin^T 



de dT 



(47) 



Finally, the particle velocity in the z direction is given by 
dz/di = vqz mil- Combining gSj-lgTI) we find for the 
classical density: 



Pair) 



2iulJo 



f{0'aA'a) 



vaZo 



Ta sin Ta 



COS 1 ( 

sin" T„ 



-I- 



(48) 

This expression diverges at the caustics where de/dT 
vanishes. It is also singular on the focal line segments 
linking cusps of caustic surfaces fSection llll C|) a.sT = kn 
holds there. 
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B. The Maslov index 

It remains to determine the Maslov index /ia(f) of 
the trajectories. This integer tracks the sign changes 
of the determinant in Eq. H44|l along the classical path 
HE 1^ US — in other words, the number of points 
along the trajectory a where the density Pa (|48l) diverges. 
In the parallel field problem, this happens (i) when the 
trajectory traverses a focal line segment {t = jn, where 
< t < r is the time of flight along the path); and (ii) 
when the trajectory is reflected at a caustic surface, i.e., 
the function ^ becomes zero along the path. Inserting 
the equation of motion ®, {Tj) into Eq. H28|) we obtain 
the equivalent condition: 



2 1 
sin^ e' cott+ - cos e' + - 



V 



t 



cos 



9' = 



(49) 



Within each cyclotron period {k — l)7r < t < kir, the 
left-hand side of this expression is a strictly monotonic 
function that approaches ±oo at the interval boundaries: 
For a given emission angle 6' there is a single solution 
to Eq. 149(1 in each period, so every trajectory eventually 
touches every caustic surface exactly once. (Note that the 
implicit equation 149|) is equivalent to the parametriza- 
tion of cos6'^^(i) H31|l obtained earlier.) 

Figure[5]illustrates that for the pair of trajectories with 
times of flight (fc — 1)t: < t < kir, ^ is negative for the 

early solution t'^'' but positive for the late solution Ty'^ 
of Eq. H21|) . The change of sign implies that the late 
path, unlike the early path, has been reflected off the 
fcth caustic surface. Since both trajectories have crossed 
[k ~ 1) focal lines and turned at (fc — 1) caustics before, 
we find for their respective Maslov index: 



(50) 



C. Caustics and a uniform approximation 

The semiclassical approximation ()48f) generally per- 
forms remarkably well, but it diverges for those trajec- 
tories that cross over each other at a caustic. To correct 
the divergent behavior, we may replace the semiclassi- 
cal contribution of the offending classical paths in the 
sum (|39|) by uniform approximations, smooth functions 
that describe the refracting electron wave near the caus- 
tic. Catastrophe Theory l4ll offers a mathemat- 
ical formalism to classify the points of the caustic set 
in terms of a small group of generic cases, the "catas- 
trophes." Several examples of higher-order catastrophes 
occur in the parallel field problem, and we list them in 
Table Inl They are located on the symmetry axis {p — 0). 

In contrast, the caustic surfaces off the symmetry 
axis {p > 0) all represent the most basic catastrophe, 
the fold. Electron waves undergoing reflection at the 
fold-type caustic surfaces locally resemble Airy functions 



Type 


Generic Catastrophe 


Caustic Surface 


tola 


Focal Line 




Cylindrical Cusp 


cusp 


Conical Lid 


hyperbolic umbilic 


Double Cusp 


hyperbolic umbilic 



TABLE IT: Catastrophe classification of caustic points in par- 
allel fields. All points of the caustic suface with p > belong 
to the simplest class, the fold; a sequence of higher catastro- 
phes is found on the symmetry axis. (Double cusp catastro- 
phes occur only in the limit E and are not discussed in 
this paper.) 



|37ll43ll44j . This suggests an approximation scheme that 
works uniformly for all p > 0: In each cyclotron period 
(k — l)7r < T < kn we replace the semiclassical contri- 
butions V'Q(r) = '0<,>(r) (|S7|l due to the two trajectories 

with time of flight T^^ in this interval (see Fig.|2) 
by a linear combination of Airy functions: 

V'itlW = lk{v) Ai(u,.(r)) + 4(r) Ai'(u,.(r)) . (51) 

To obtain the uniform approximation we then sum 

V-'ink'") '^^'^^ ^ cyclotron orbits k. 

We choose the auxiliary functions in Eq. (|51|l so that 
after asymptotic expansion of the Airy functions i the 
leading term matches the combined semiclassical contri- 
butions of the two classical paths in the kth cyclotron 
period: The argument of the Airy functions Wfc(r) only 
depends on the difference in the classical action of the 
two paths, AWfc = 'Wk,> - Wfc,<: 



Mfe(r) 



(52) 



while the "amplitudes" 7fe(r) and Sk{v) are functions of 
the classical densities ^^(r) (|48ll along these trajectories: 



7, = (-l)'=V^|z,,|V4eiw.A-W4(^^+^;^) ^ 



(53) 



(54) 



Here, Wk ~ \{y^k,> + W'fc,<) denotes the average action 
for both paths. A similar identification scheme applies in 
the case of classically forbidden motion, where we use a 
tunnehng trajectory to fix Wfe (r ) , 7^ (r ) , and (5^ (r ) instead. 
These auxiliary functions, and therefore the uniform ap- 
proximation (|51|l . behave smoothly at the caustic. 

As noted before, the approximation is only valid for 
fold catastrophes, and, like the semiclassical solution, 
fails near the focal lines, cusps, and umbilic catastrophes 
that occur on the symmetry axis. These regions require 
a separate treatment. 
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V. THE QUANTUM SOLUTION 

In the quantum description, the motion of the electrons 
in the external potential A(f), $(f) ^ is controlled by 
the Hamiltonian H : 

H(f , p) = _ . p - dv, p) = — (p - qA(f ))' + qm . 

(55) 

As in Sectional we again introduce dimensionless vari- 
ables. A natural energy scale e for the perpendicular 
motion is given by: 

e = E/ihwL) , (56) 
and Ti. takes on a particularly simple form if we set: 



n V h 



y = V2ey, (57) 



For the motion in the field direction, we define an inverse 
energy scale /3: 



f3 



2m si/3 



and a dimensionless coordinate ^: 



C = PqSz = 2 — 
r] 



(58) 



(59) 



Because the motion along the symmetry axis and in 
the perpendicular x-y plane are independent, the Hamil- 
tonian H 1)55(1 is separable: 



(60) 



and the partial Hamiltonians commute, [7Yj^,7i||] = 0. 
Using dimensionless variables, they read: 



1 



i(C9. 



2 ' v'] 



and 



H||(£,p.) = --[9^2 + C] 



(61) 



(62) 



A. Electronic eigenfunctions 

The Hamiltonian ni_ (EU governs the cyclotron mo- 
tion of an electron confined to the x-y plane in a uniform 
magnetic field. It differs from the Hamiltonian of a two- 
dimensional harmonic oscillator only by the middle term 
— wlL^ which is constant for radially symmetric eigen- 
states with fixed angular momentum quantum number 
II. The stationary states are thus oscillator eigenstates 
that are shifted in energy; we express them in terms of 




-3-2-1 1 2 

Angular Momentum [in units of h] 




e -o 



Lz<0 



Lz = 



Lz>0 



FIG. 6: Top: The spectrum of Ti.±. The eigenstates can 
be ordered by their angular momentum Lz — ^ih, and their 
energy Ei, = {2u + l)hLUL- Their distribution discretizes the 
classically allowed range -B > 0, -E > —2ujlLz (red area). 
— Bottom: Corresponding classical orbits in the x-y plane. 
The energy E — 2mLo\R^ determines their radius R, while 
the angular momentum Lz contains information about their 
position relative to the origin: The orbit encircles the origin 
for Lz < 0, while the origin is outside the orbit for Lz > 0. 
Lz vanishes for orbits passing through the origin (like the 
trajectories of Section^. (The magnetic field points toward 
the reader.) 



Laguerre polynomials i^J(pf) |4J|, where = + is 



the scaled radial distance 



tt/i {n -)- 1^1)! 



pl^lLlfl(p2)e-^'/2eiM0 



(63) 

For /X > 0, but not for negative /x, their eigenenergy 
depends only on the principal quantum number n — 
0,1,2,...: 



(2n-|- \)hwi 



(m>o) 



{2n^2\ix\^\)hjjL (Ai<0) 



(64) 



The spectrum is discrete, and the states occupy equally 
spaced, infinitely degenerate Landau levels = {2v + 
l)hujL 1131 . The unusual skewed spectrum can be ex- 
plained in terms of classical cyclotron orbits (see Fig.lHJ. 



The Schrodinger equation for motion along the z axis, 
— '^||I21|(C; -E'li) = 0, is essentially Airy's differential 
equation |47| , and the eigenfunctions can be expressed as 
regular Airy functions |44| : 



iikii(C;- 

|47l |. an 

■ fun< 

xLg(C;£^ii) = Ai[-(/?i?||+c)] 



(65) 
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For our purposes, we will also need a Hankel-type solu- 
tion of the differential equation that behaves like an out- 
going wave in the asymptotic limit z ^ oo 0, |23, E3| '■ 



contribute to the series (|69|l (all other eigenfunctions 
XuAp'A') vanish for p' = 0). Because L?(0) = 1 [3, 
we find the simplified expansion jl^ |33, : 



xL(C;i?ii) = Ci[-(/3i?||+c)] 



(66) 



Here, Ci{u) — Bi(7i) + i Ai(u) is a linear superposition of 
the standard Airy functions 44] . 



B. Construction of the Green function 

The wave function 'ijj{r) generated by a stationary, 
pointlike, isotropic source of electrons with energy E at 
position f ' is a solution of the inhomogeneous Schrodinger 
equation p^ : 



[E - n] V'(f , r';E) = CS^^^ (f - f') 



(67) 



The parameter C describes the strength of the electron 
source. The solution ip of Eq. H67I) is then a multiple of 
the energy Green function G{r,v']E) which is formally 
the resolvent operator (i.e., the inverse oi E — Ti) in po- 
sition space p8| : 



G{v,r'-E) 



lim 



(68) 



The infinitesimal parameter fj ensures that the Green 
function asymptotically shows outgoing-wave behavior. 

The quantum mechanical solution therefore requires 
knowledge of the Green function. In three-dimensional 
space, explicit expressions for G'(f , f ' ; E ) a re available 
only for a handful of simple potentials [13, |43, ISfl 1^ > not 
including the parallel fields Hamiltonian Ti (I6UII . How- 
ever, G'(f , f'; E) can be written as a series of readily avail- 
able one-dimensional Green functions [HJI: We first ex- 
pand Eq. into the complete set of eigenstates Xn.^ 
()63ll . (|64|) of the "magnetic" Hamiltonian 7ij_ = ?i — 7i|| : 



G{v,v'-E)=Y,xiAp,<f)xi^{P'Ar 



X lim ( z fi? - 



(69) 



Comparison with Eq. (|68|l shows that the second line is 
actually the outgoing- wave Green function for an electron 
with effective energy E\^fj, ^ E — E^ moving in the 
one-dimensional linear potential 7i|| H62f) . To obtain it 
in explicit form we match the regular and outgoing-wave 
solutions for a uniformly accelerated electron H65|l . H66fl 
at the source z = z' [T^ : 



G|| (C, C'; E\\)= -^P'q£ Ai(-u<) Ci(-u>) 



(70) 



Here, u< ~ (3E\\ +mm{(X') and u> = -|-max(C,C')- 
We now assume that the source is located at the origin, 
f' = 0. In this case, only magnetic eigenstates with ^ = 



G{v;E) 



mB 



X M[-I3{E - En)] Ci[-[3{E - E^) - C] , (71) 

where En = (2n + l)hujL is the energy bound in the 
cyclotron motion (|64|l . The result holds for z > 0; for 
z < 0, the arguments of the Airy functions have to be 
exchanged. 

For large n, the energy = E — En available for mo- 
tion along the z axis is negative, the electrons are forced 
to tunnel from z' = to z, and the corresponding term in 
the expansion (|71|l is exponentially suppressed. There- 
fore, this series generally converges rapidly, and permits 
accurate determination of the quantum solution G{r; E). 
(This even includes tunneling sources with E < Q that 
have no classical counterpart. The only exception is the 
plane z = where the series does not converge absolutely. 
This is not surprising because G(f ; E) itself diverges as 
f approaches the origin, G{r;E) ^ —m/{2nh'^f).) 



C. Currents generated by the source 

The stationary Schrodinger equation (|67|l contains 
an additional inhomogeneous term, the "point source" 
C(5'-^-'(f). Defining the current density distribution in the 
scattering wave ip{r]E) in the usual manner, 

j{r;E) = -lm[rV^]-^^\^\\ (72) 

we integrate j(f; E) over a surface dfl enclosing the origin 
and find the total current J{E) emitted by the source: 



J{E)=(f j(f;£:) •dn(f) = / d-VV-j(f;£;) 
Jan Jn 



(73) 



Solutions of the usual homogeneous Schrodinger equation 
fulfill V • j = and cannot carry a net flux through any 
closed surface. The presence of the source modifies the 
equation of continuity, however |27l |: 

V-j{r;E) = -| Im[C*<5(3) (f . e)] , (74) 

and because ^p{r;E) = CG{r;E), the total current is 
proportional to the imaginary part of the Green function 
at the source location 



JiE) 



|Cp liniIm[G(f;£;)] 



(75) 



The parallel field Green function H71|) yields a series 
representation for J{E). We obtain with L„(0) = 1 [i^ : 



JiE) 



2mB 



\c\'j2H-i^iE-En)y 



(76) 



n=0 
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J{E) is experimentally accessible as the total s-wave 
photodetachment cross section. Similar results for the p- 
wave cross section have been found using Fermi's Golden 
Rule [Ulilil. 



VI. EXAMPLES 

To illustrate the results of the two previous sections, 
we now plot and discuss radial profiles of the integrated 
density n{p) = 27rp|'!/'(/5)P for three cases that sample 
different regimes of motion in parallel fields. (We use 
integrated densities because the electrons tend to cluster 
around the symmetry axis (see Fig. indeed, for p — > 
the semiclassical approximation nsdp) approaches a 
constant if a focal line is present.) All results shown are 
normalized to a total flux of one electron per second. 



A. Moderate Fields 

We first return to the model problem used in the in- 
troduction {S = 15 V/m, B = 0.02 T, E = 10"" eV, or 
77 = 7.908, e = 86.38; see Fig.P) as an example with mod- 
erate parameters — semiclassically, at most eight trajec- 
tories contribute, while the quantum mechanical solution 
converges after summation over about fifty Landau lev- 
els. Figure [7| shows details of the radial density profile 
n(/5) at a detector distance z — SO fim. The semiclassical 
result nsc(p) obtained including complex "tunnel- 
ing trajectories" (green) is an excellent approximation 
to the exact solution n{p) (blue) given by the quantum 
Green function H71|l . except near the four intersections 
with caustic surfaces (red dotted lines) . The uniform ex- 
pansion (|51|l into Airy functions (not shown) is virtually 
indistinguishable from the exact result. Its typical rela- 
tive error is much smaller than 1 % over most of the radial 
range, except in the immediate vicinity of the symmetry 
axis: Near a focal line, the densities Psc{p) and Puni(/5) 
diverge like p~^. Therefore, nsc{p) and nuni(/5) approach 
a constant value as p — > 0, while n(p) vanishes linearly. 
A proper treatment would require a different uniform ap- 
proximation valid near focal line singularities. 

Another feature of Fig. [7| is the "metastructure" vis- 
ible in the interference pattern as a modulation in n{p) 
around p 2 pm. Its length scale (~ 0.3 /im) falls in 
between the period of individual oscillations in the den- 
sity (~ 0.08 /Ltm), and the range of the classical caustic 
structure (~ l/xm). An approximate symmetry of n{p) 
with respect to the location of the snake- type closed orbit 
(which intersects the detector plane at pco = 1.945 ^m) 
is conspicuous. We already pointed out the role of closed 
orbits in organizing the spatial interference pattern in the 
electronic density (see Fig. ,6^J . 




Radial Distance ^ [fim] 



FIG. 7: Integrated density profile for a point source of elec- 
trons (normalized flux 1 s~ , energy E = 10"" eV) in paraUel 
fields £ = 15 V/m, B = 0.02 T (r? = 7.908, e = 86.38), at a 
distance z — 30 pm from the source (cf. Fig. 0. Blue curve: 
Quantum result 17111 . green curve (dashed): Semiclassical ap- 
proximation 1391 , including tunneling trajectories. The WKB 
approximation diverges at the caustics (red dotted lines) . The 
red arrow marks the position of the "snake" closed orbit. 



B. Weak Magnetic Field 

The addition of even a weak magnetic field to an elec- 
tric field fundamentally alters the structure of the tra- 
jectories and the resulting interference pattern. A pure 
electric field produces a "fireworks" structure of trajecto- 
ries, and just two paths arrive at a point on the detector 
[itI l2^ . In contrast, with any magnetic field present, 
the emitted electrons initially spread outward, but the 
magnetic field eventually turns them around and rcfo- 
cuses them onto the symmetry axis, where they form a 
focal line segment that joins opposite cusps of adjacent 
caustic surfaces f Section llllll . For small ry, the caustics 
barely overlap and form a sequence of narrow, diamond- 
shaped focal regions with four-path interference. They 
are centered around Zk — (kn)'^ /i] and have an approx- 
imate radial width E/{q£). Outside the focal regions, 
two-path interference prevails; the behavior matches the 
purely electric case, except that the cyclotron motion 
causes the radius of the pattern to oscillate. With in- 
creasing distance z, the radial electronic distribution ex- 
pands and contracts, and the current is forced through a 
sequence of hourglass-shaped constrictions (see Fig. |3Ia) 
and Ref. jH). 

Figure |S1 depicts radial density profiles at maximum 
extension (top) and at the first "bottleneck" (bottom) 
for the set of parameters £ — 116 V/m, B = 1 mT, and 
E = 10-* eV, or 77 = 0.051 and e = 1728. In the first case 
{z — 3.255 mm), we obtain a simple two-path interference 
pattern. The semiclassical approximation (green) and 
the quantum solution (blue) agree very well except near 
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0.2 0.4 0.6 0.8 

Radial Distance ^ [|j,m] 



FIG. 8: Integrated density profiles for a normalized point 
source in fields f = 116 V/m and S = 1 mT (electron en- 
ergy: E = lO""* eV, r) = 0.051, e = 1728). Top: Distribu- 
tion close to the maximum radial extension (detector distance 
z = 3.255 mm), bottom: Distribution at the "bottleneck" 
constriction [z = 13.02 mm; note the change in scale for the 
p ax;is). The uniform approximation 15111 (red) is virtually 
indistinguishable from the quantum result II71II (blue) , except 
near p = Q (see detail). The semiclassical approximation in- 
cluding tunneling paths 1391 (dashed green curve) diverges at 
the caustics (red dotted lines). 



the outer turning point at p = 67 A /^m. Because no focal 
line is present, the WKB approximation holds even for 
p — 0. At the center of the constriction (z = 13.02 mm), 
a more involved oscillation pattern emerges as four paths 
interfere with each other in the sector of classically al- 
lowed motion. The semiclassical approximation fails near 
the double caustic at 0.862 pm, and also close to the 
symmetry axis due to the presence of the focal line, but 
performs well otherwise. To repair the divergence at the 
caustics we resort to the uniform approximation (red) 
which is indistinguishable from the quantum result ex- 
cept near p = 0. Note that the oscillations in the density 



profile extend beyond the turning point; they indicate 
interference of the two tunneling trajectories. 

For weak magnetic fields, a large number of terms 
(~ 10'^) must be summed up to achieve convergence in 
the series for G{r;E) H71|) . (Note that they cancel each 
other almost perfectly everywhere outside the classically 
allowed region; this cancellation over a large region is 
especially remarkable near the bottleneck.) In contrast, 
the semiclassical method relies on at most four trajecto- 
ries. Here, the uniform approximation (|51|I - H54|I is a very 
efficient scheme to evaluate the wave function. 



C. Strong Magnetic Field 

Finally, we take a look at the opposite limit of a strong 
magnetic field B. Figure displays the density profile at 
a distance z = 30 ^m for £ = 15 V/m, B = 0.5 T, and 
E — 10~^ eV, corresponding to dimensionless variables 
T] = 197.7 and e = 3.455. The situation is most natu- 
rally described in terms of the quantum result n{p) (blue 
curve) — effectively, only two open "channels" contribute 
to the series for G{r';E) H71|l . In contrast, up to 128 
classical trajectories connect the source to a point on the 
detector. The 64 intersections with caustics (marked by 
red bars) are spread over the entire range of classical mo- 
tion, and they are most densely clustered near its outer 
limit p = 1 (or p = 0.135 ^m). These caustics render 
the primitive semiclassical result (|39|) virtually worthless. 
However, it is striking that the uniform semiclassical ap- 
proximation niini(p) still performs well over most of 
the range of p (except near p = 0, where it fails due to 
the presence of focal lines). This is an example where we 
may reasonably expect this approximation to fail, but it 
gives a surprisingly good result. We conclude that the 
expansion into Airy functions (dashed red curve) is a re- 
liable approximation over a wide range of parameters. 



VII. CONCLUSION AND OUTLOOK 

We studied the dynamics of monoenergetic electrons 
emitted from a pointlike, isotropic source placed in ho- 
mogeneous, parallel electric and magnetic fields £, B. 
The problem is simple enough to permit a full quantum 
solution in the form of a series for the Green function, 
Eq. (|71|l : likewise, analytic expressions are available for 
the classical electron motion, and we used them to de- 
velop semiclassical and uniform Airy approximations to 
the wave function. 

By tuning the two "control parameters" rj ~ vqB/£ 
and e = E/(fiWL) that characterize the system in the 
classical and quantum domain, both the number of in- 
terfering classical trajectories and the number of open 
quantum scattering channels can be adjusted indepen- 
dently. This is accomplished by simply varying the elec- 
tron energy and/or the external fields. 



19 




Radial Distance p [|im] 



FIG. 9: Integrated density profile (normalized flux) for £ = 
15 V/m, B = 0.5 T, and E = IQ-* eV {rj = 197.7, e = 3.455), 
at a detector distance z — 30 ^m. Blue curve (bold); Quan- 
tum result 17H . red curve (dashed): Uniform approximation 
1511 1. green line (thin): WKB approximation 1391 with tun- 
neling contributions. The sequence of bars (red) marks the 
position of the 64 caustics. The purple curve indicates the 
number of interfering classical trajectories N{p). 



We verified that the patterns seen in the quantum 
density distribution are naturaUy explained in terms of 
the caustic surfaces and closed orbits of the correspond- 
ing classical motion. Using a temporal parametrization 
technique, we examined the set of caustic surfaces and 
found a surprisingly complex classification scheme fea- 
turing seven distinct types of caustics. With its help, we 
learned how the caustic pattern systematically evolves 
with increasing rj. 

The present state of our study suggests future work in 
two directions: While we achieved a complete description 
of the caustic structure for p > 0, we did not examine in 
detail the wave functions near the focal lines, cusps, and 
umbilical points on the p = axis. In the future we plan 
a systematic study of these catastrophes and their proper 
uniform approximations in the parallel field problem. 

Pointlike sources of electrons are almost perfectly re- 
alized in near-threshold photodetachment of negative 
atomic ions. Photodetachment experiments in parallel 
fields should therefore provide practical access to a sys- 
tem where the interference of matter waves can be pre- 
cisely tailored. For the field strengths typically used in 
photodetachment microscopy, the interference patterns 
tend to be very small, however (with diameters in the 
/im range), which makes them hard to observe. To over- 
come this difficulty, we will examine non-uniform field 
configurations that "flare out" from the source and effec- 
tively magnify the density distribution on the detector. 
Given the excellent accuracy of the semiclassical and uni- 
form approximations reported in this paper, we expect 
that semiclassical models likewise yield good predictions 
for the interference patterns in these more sophisticated 
field geometries. 
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APPENDIX A: TEMPORAL 
PARAMETRIZATION OF CAUSTIC SURFACES 

In Section UTIl we described the caustic structure 
emerging in parallel fields, and its evolution with 
changing value of the parameter rj. The results pre- 
sented there were based on the temporal parametrization 
[/9±(i), 2:±(t)] of the caustic surfaces. However, for 
the sake of clarity, we furnished no proof in the main 
text. We now discuss the mathematical properties of the 
parametrization and show how they may be exploited 
for a complete characterization of the caustics. While 
our argument employs only elementary mathematics, the 
complete description is surprisingly complex. 

1. Proof of the algorithm in Section llll Bl 

We first give a proof for the algorithm shown in Fig- 
ure 13 The existence of solutions in H29|l is contingent 
upon the conditions laid out in Eqs. ()32a|l . (j32b|l . We 
consider only t > 0. 

a. T < 

We start with the case r = tant < 0. Then, A{T,t) > 
> (|30() is always positive, and we can restrict our 
attention to the second inequality. 

Since t + \fA > t + |t|, the second condition in (|32b(l 
is always satisfied for the (-I-) solution. We show next 
that Eq. (|32b|) is fulfilled for the (-) solution only if 
t < I holds. For the proof, we compare both sides in 
this inequality: 

(t~Tf <{t-^fA)' . (Al) 
When multiplied out, this is equivalent to: 

2 

2t^^A<{2t-^^T + =B{T,t) . (A2) 

To facilitate comparison, we first square both sides and 
form their difference: 

A(t, t) = B^- 4fA = if {if - 4f) (l - ' , (A3) 

from which it follows (barring t — t) that sgnA(T, i) — 
sgn(| — t). The condition HA2II is then equivalent to 
A > and B > 0. In the interval t > ^, the first 
inequality is never satisfied, and hence the (— ) solution 
is complex there. For t < ^ and r < 0, however, both 
inequalities hold, and the (— ) solution is real. 
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FIG. 10: The auxiliary lunctions A(T,t) ( ereen, solid 

curve) and B{r,t) iA2l (blue, dotted line) for t < f. In 
the red intervals, the parametrization H29^ describes a pair of 
physical solutions. 



b. r > 

Now we assume t > 0. Using the same line of argu- 
ment, we show that for the (— ) solution, Eq. (|32a|l implies 
2tVA > B, or equivalently, A < or B < 0. Therefore, 
the condition is fulfilled if < > |. Since A{T,t) > al- 
ways holds in this regime, we find that the (— ) solution 
is real for t > f and tani > 0. For the {+) solution, 
we find the condition —2t\/A > B, which is identical to 
the inequality pair A > and B < 0. For t > |, they 
are never simultaneously fulfilled, so the (+) solution is 
complex for T > [fi^j ]. 

However, these inequalities show that for r > 0, both 
solutions will be real in the interval t < | if A{t, t) > 
and B{T,t) < 0, and otherwise complex. A{T,t) is a 
quadratic function of t that formally has roots: 



c. t = rj/2 

We finally examine the behavior at the boundary t — 
^. Above, we inferred that the solution space is degen- 
erate if either (r = tan ^ < 0), or {A{T,t) > and 
B{t, t) < 0). Because 

A{T,t) = (r;-tan|)' , B(t, i) = 77(77- tan |) , (A6) 

the latter condition holds for tan § > '7- the intervals 
< tan § < 77, only the (— ) solution is real. 

2. Irregular points on caustics 

A point on the caustic (^o, ^o? Po) is regular if the func- 
tions /3±(i), z±(t), and either p{z) or z[p) have continuous 
derivatives (C°°) at this point. Otherwise the point is ir- 
regular. The caustic surfaces in space must be smooth 
unless one or both of the functions p± (t) , z± (t) become 
irregular, or both derivatives dp±{t)/dt and dz±{t)/dt 
simultaneously vanish. 

The parametrization (|29|) has irregular points at (a), 
(b) to = {k — ^)Tr so T — tani diverges; (c) to = ^ so 
the radicand in p±{t) vanishes for either the or the 
(— ) solution; (d) the minimal values to = for which 
A{ta.iito,to) — and B{tanto,to) < (see above); 
(e) to = kn so T = (see Section llIIC|) . These irregular 
points occur for general values of 77, and each case occurs 
at most once in every cyclotron period (fc — l)7r < t < 
kn. For isolated values of rj, two of these conditions may 
coincide. We will address this issue separately below. 

The behavior of the caustic surface and its temporal 
parametrization near the irregular points (to, zo, po) may 
be inferred from a local expansion of Eq. 129|l into powers 
of the deviations St — t — to, Sp = p — po, Sz — z — zo [63| . 
The results of this tedious yet straightforward calculation 
are assembled in Table IIIII Here, we proceed with an 
interpretation of our findings. 



Within the interval r< < r < t>, A{T,t) is negative. 
Similarly, B{T,t) is a linear function of r with negative 
slope; B{Q,t) — rf, and B{T,t) remains positive for < 
T < tb = trf' /{rf' — 2t^). In both intervals, all solutions 
are complex. At the zero tb of B(t, t), however, A{tb, t) 
is negative: 

/4t^ \ 

A{TB,t) = rl[—-l) . (A5) 

Therefore, the conditions A(t, <) > and B{T,t) < arc 
simultaneously satisfied when, and only when, t > Ty 
(see Figure [TU|l . 



a. Maxima of p: to — {k — ^)tv > 2 

This is the simplest case. Here, r — tant diverges, 
and as t increases through this point we switch from the 
(— ) solution to the (-I-) solution (see algorithm in Fig.|21l. 
The caustic is a smooth curve in (p, z) space. This point 
represents the maximum radial extent of the caustics. 

6. Maxima of p vs. endpoints of closed orbits: 
to^{k- i)7r < f 

Here, the caustics are smooth curves in (p, z) space but 
the parametrization is tricky. Both (-I-) and (— ) solutions 
exist, but are discontinuous: Passing through the diver- 
gence of tant, the (— ) curve [p_ (t), z^ (t)] joins smoothly 
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Local Behavior of Caustics 
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to = [k- i)7r 


to 


1 


sz^'^^^^:+'^ St 
277 

''--2tJW+^)'' 


Maximum radial extension of caustics; sin- 
gularity in the parametrization only; switch 
from (+) solution tor ot > to (—) solution 
for St<0 


to = (fc- |)7r 
77 > (2fc - l)7r 


£0 


77 V ' " 


,__to[.^+4(l-t.1)] 
277 

5p ~ 5z 


Singularity in the parametrization, existing 
only for to < -1; switch from (— ) solution for 
ot > to (+) solution tor Ot < 0. Corre- 
sponds to a "snake" -type closed orbit 


to = i 
77 7^ fcyr, 

r] + Vcrit 
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Sz ' ^ St 

77 — tan ^ 

X 2 4 ■ 2 V s 
op ~ — sm — oz 


Uphill turning point of parabolic shape; sin- 
gularity of the parametrization only. St is 

nnt^itivp for Oik — 1 \it <^ T) ^ ■ anH npfffi- 

five for 77^111 < 77 < 2A;7r 


J? > Vcvit 


tan to 


\/ sin^ to — 772 cos^ to 


tan to at 

to 

77sin(2to) , 

op ~ Oz 

2poto 


Minimum time of flight along caustic; sin- 
gularity of the parametrization only. Switch 

ilUlli \ \ ) bUiLlLiUli iUi uZ ^ U LU I 1 OUIULIUII 

for Sz>Q 


to = fcvr 


kTv{kn + 77) 
V 





3(77 + 2fc.) 
277 

2 

27fc7r(77 + 2fc7r)2 


Cylindrical cusp, opening towards the source 
(Sz < 0). Physical singularity of the 
solution. St must be negative 


to = fcvr 
77 7^ 2kn 


kTv[kn — 77) 





Sz^ --^i^ St 
277 

.2 8772 3 
^ 27fc7r(77 - 2fc7r)2 


Cylindrical cusp, opening in field direction 
(Sz > 0). Physical singularity of the (— ) 
solution. Sign of St depends on to: St > 
for to > f , <5t < for to < f 




Irregular Points for Special Values of 77 


to 


zo 


Po 


Local Behavior of Caustics 


Remarks 


to = kn 
77 = 2kn 


V 
4 





Szr-. 1- St'' 

kn 
Sp ~ Sz 


Umbilical point in the (— ) solution; physical 
singularity of the caustic set; two surfaces 
(distinguished by the sign of St) form a dou- 
ble cylindrical cone with opening angle j 


to = i 


4 





Sz'' ~ 77<5t 

>c 2 4772 
1 + 772 


Bifurcation point in the temporal 
parametrization (caustic remains parabolic); 
solution at top of caustic switches to (-I-) 
branch; extreme sensitivity to changes in 
St>0 


to = (fc — |)7r 
77 = (2fc - l)7r 


V 
4 





(52 ~ -St 

2 4 
5p ~ — 5^ 

V 


Creation of a "snake" -type closed-orbit; sin- 
gularity of the parametrization only. At the 
top of caustic (which remains parabolic), so- 
lution switches to (— ) branch 



TABLE III: Irregular points of the temporal parametrization [z±{t), p±{t)^ Ij29|l . together with local leading-order expansions 
for the caustic surfaces. (The minimal times t\^^^ and critical parameter values 77^^;^ are defined in Eqs. Ij.'M^ and iPSJO 



to the curve [p+(t), , and vice versa. As St in- 

creases continuously along the caustic, if we switch from 
the (— ) solution to the solution, the caustic has a 
radial maximum, just as in case (a). If we switch from 
the (-I-) solution to the (— ) solution, the caustic curve 
is again smooth, and the irregular point represents the 
endpoint of a "snake" -type closed orbit where p = z = 0. 



c. Top of caustic: to — ^ , rj ^ kir and 77 7^ 77*^ 



This is the topmost point of the caustics (zq = ^4)) 
where z{p) has a quadratic maximum. The parametriza- 
tion is again tricky: dz is proportional to St, and Sp 
is proportional to |(St|^/^. This irregular point in the 
parametrization occurs in the (— ) solution if tan ^ < 77, 
and in the branch otherwise. Further analysis (see 
Table shows that St must be positive if 2(fc — l)7r < 
?7 < r?ci.it: and negative if 77^^;^ < 77 < 2kn. The critical 
values of 77 are given by rj^^.^^ = tan(?7^^;j/2) 
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while a more complex law connects 5p and 5t: 
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FIG. 11: Cusps and cones in caustic surfaces near the irreg- 
ular point (po = 0,20 = — f). In tiie umbilic case to = tt 
and 77 = 27r, a pair of surfaces asymptotically coalesces into 
a double cone of opening angle 45° ("early" caustic 5t < 0, 
green line; "late" caustic St > 0, blue curve). In contrast, 
for to = 2n the generic value — 8n/3 leads to a cylindrical 
cusp (red dotted curve, St > 0). A log- log plot of Sz ver- 
sus Sp (inset) illustrates the different power laws 5p ~ Sz and 
Sp^ ~ C{r])Sz^ guiding these surfaces, respectively. 



Sp: 



7]_±2fc7r^^3 
rjkTT 



(A7) 



For a meaningful result, the sign of St in Eq. I|A7|I must 
be properly chosen: We always have St < in the 
solution; in the (— ) solution, St is positive (nega- 
tive) for to > I {to < §). (The special case r] — 2k-K 
will be discussed later.) Eliminating St, we find that 
the caustic surface is locally parametrized in the form 
5p\ ± C{r])5z\ = where C{ri) > is a positive coeffi- 
cient. This is the characteristic functional form of a cusp 
singularity js^. The caustic surface follows from rota- 
tion of the cusp profile around the symmetry axis and 
has been called a cylindrical cusp 

The (-I-) and (— ) solutions produce cusps that point in 
opposite directions: The cusp in the (+) solution points 
downhill while that in the (— ) solution points uphill in 
Fig. ^ Between each pair of cusps at z^(fc7r) a focal 
line segment extends along the z-axis, where trajectories 
from all values of the polar angle (f)' are focused onto 
p = 0. These trajectories simultaneously arrive at the 
z-axis after k complete cyclotron orbits. Their emission 
angles 9' vary between (downhill) for and tt (uphill) 
for Zq. Depending on k and 77, the cusps at Zq and z^ 
may or may not belong to the same caustic surface. 



d. Minimal times: yl(tanto,io) = 

Again, the caustic curve is smooth but the 
parametrization is irregular. If t < ^ and rj > 77^^.;^ then 
in each cyclotron period (fc — l)7r < t < kir there are 
two times at which yl(tant, t) = 0, but in the vicinity of 
the smaller root, no physical solution exists for other rea- 
sons (see previous section). Therefore, only the larger of 
the two roots, the minimal time t^^}^ l|34|) . is an irregular 

point. As t increases through t^^/^, the radicand A(t, t) 
(|30|l passes from negative to positive, and a pair of so- 
lutions [z±(t), p± (i)] springs into existence. The time to 
hence marks a local minimum of the time of flight from 
the source to the caustic. From this point, the caustic is 
described by the (— ) solution for positive Sz, and by the 
(-I-) solution for Sz < 0. In the vicinity of the singularity, 
small variations St > oi the time of flight cause large 
shifts in position. 



e. Cusps: to — kn, rj 7^ 2fc7r 

In all previous cases, the irregularity occurs in the time 
parametrization but the caustic itself is a smooth curve 
in (p, z) space. This case is different. When to = kir, 
T vanishes, and the irregular points lie on the z-axis at 
= ^(k'!T±ri). As stated in Table lllll in the vicinity of 
these points the relationship between Sz and St is linear. 



3. Special irregular points of caustics 

The above irregular points are "generic" in the sense 
that they occur at almost all values of rj. Three extra 
cases arise only for special values of 77, when two generic 
irregular points coalesce. 

a. Cylindrical hyperbolic umbilic: to = f ~ 

If ry is an integer multiple of 27r, then the time to ar- 
rive at the top of the caustic, to = 2 2c), coin- 
cides with a completed cyclotron orbit (case 2d). As rj 
increases to 2fc7r, a cusp moves up to collide with the 
top of the caustic. The cylindrical cusp and the smooth 
dome are changed into a double cylindrical cone having 
two branches, both meeting the z-axis at an angle of j. 
These branches are both described by the (— ) solution 
which is continuous at to = kn; they differ in the sign of St 
(note that Sp ^ Sz St^ /kir). Figure ITTI shows the caus- 
tic for T] = 27r. In Catastrophe Theory jH EE El El , 
this represents the generic catastrophe in three control 
parameters {p,z,ri) known as the hyperbolic umbilic. 

b. Bifurcation points of the parametrization: 
to = %,r) = rj^^it 

This irregular point occurs at the critical values of 77 
given by Eq. Ij^Sfl when cases (2c) and (2d) coincide. As 
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p [dimensionless] 

FIG. 12: Bifurcation and creation of a "snake"-type closed 
orbit in the first caustic surface {k — 1). The figure displays 
a family of graphs of the function t{p) with increasing values 
of ri (bottom to top): For rj < T^^^riti *(p) has its minimum 
on the symmetry axis at the top of the caustic, t{0) = § 
(case 2c), and p{t) grows monotonically in < p < 1 (bottom 
curves, = 2, 2.17, corresponding to the late variant [SI] 
of the "teardrop" caustic. See also Fig. |1]) The behavior 
changes at the critical value T^J^t ~ 2.33112... (bottom heavy 
curve): For i] > r/criti p = is a local maximum of t{p), while 
the minimum t = t^^l^ (case 2d) is now attained off-axis. We 
are at the intermediate-time variant [Si] of the "teardrop" 
{ri = 2.50, 2.67, 2.83, 3). The extrema of t{p) are located 
on the dashed green curve ("pitchfork bifurcation"). — At 
r] = TT (top heavy curve), the first snake orbit with time of 
fiight t = f- (dashed purple line) detaches from the symmetry 
a:xis |3, and we enter the early regime [Se] of the caustic (77 — 
3.33, 3.5, top curves). — Red and blue segments represent 
the (-I-) and (— ) solution in Eq. 12911 . respectively. Symbols 
indicate irregular points of the parametrization: Half-integer 
orbits (□, cases 2a, 2b), the top of the caustic (★, case 2c), 
the minimal time t^J^}^ (0, case 2d), the bifurcation point (v, 
case 3b), and creation of a closed orbit (A, case 3c). 



r] increases through rj^^^, keeping always to = ^, then 
according to the algorithm in Fig. O the (— ) solution 
continues, but the (-I-) solution suddenly becomes real: 
In the process, a pair of real roots of A(tanto,io) — 
is created in the fcth cyclotron period, and the larger of 

(k) 

these roots t^(^ presents the minimal time of flight along 
the caustic (case 2d) . Concurrently, the top of the caustic 
(io = 2 ' 2c) changes from a local minimum of t into 
a local maximum (see Table . 

For fc = 1, the transition is illustrated in Fig. El In- 
stead of the possibly double- valued function p± (t) , we ex- 
amine the inverted function t{p) which is locally smooth 
(we limit ourselves to the top half of the uppermost 
smooth dome caustic). We now look at the position of 
the extrema of t{p) for varying values of 1] (green dot- 
ted line) that correspond to the singularities (branching 
points) of p±{t). For ry < jy^rif ^^'^ minimum of t{p) is 
on the symmetry axis (the top of the caustic), while for 

V > ^criti the minimum occurs off-axis (at to = and 
p — is a maximum of t(p). The curve given by the irreg- 
ular points of the parametrization [to^rj), po{r])\ (green) 
has the shape characteristic for a "pitchfork" bifurcation. 



c. Creation of a closed orbit: to ~ ^ = {k — ^)it 

Here, cases (2b) and (2c) hold simultaneously. At this 
irregular point, the return time for the antiparallel closed 
orbit AT = rj — {2k — 1)tt becomes an odd multiple of 
the cyclotron period. This is exactly the return time for 
the fcth snake-type closed orbit which detaches from the 
"uphill" closed orbit at this value of 77 0. Because a 
switch from the (-1-) solution to the (— ) solution takes 
place at the endpoint of a snake closed orbit (see the 
previous discussion of case 2b), the parametrization of 
the caustic increases in complexity when 77 exceeds the 
critical value (2/c— l)7r. (See also Fig. 1121) The parabolic 
shape near the top of the caustic is not affected, however. 
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Near-threshold photodetachment into s-waves is the 
commonly observed case 28], with H being an 

important exception P]. A systematic way of dealing with 
multipole sources is developed in: C. Bracher, T. Kramer, 
and M. Kleber, Phys. Rev. A 67, 043601 (2003). 
In contrast, the charge moves uniformly with an angular 
velocity cue = 2ujl if the center of the orbit is chosen as 
the origin of the coordinate system. 
[57] The conjugate complex roots of Eq. (I21II formally corre- 
spond to "tunneling trajectories" in complex space-time. 
Because they describe evanescent waves in the classically 
forbidden sector of motion, they play an important role in 
refining the semiclassical solution near caustic surfaces. 
For a detailed analysis of the tunneling trajectories in a 
purely electric field, see Ref. [l^ . 
[58] To achieve agreement with the Green function we intro- 
duce a (physically irrelevant) overall phase factor (—1). 
We also omit contributions of order 0{h). 
For simplicity, we ignore the electron magnetic moment 
PL that merely leads to an effective shift ±^giiBB in the 
electron energy, depending on its spin orientation. 
The imaginary part of the Green function ip = 
Im[G(f ; E)] is a solution of the ordinary Schrodinger 
equation [E — TL]^' = 0. Unlike G{y; E) itself, it remains 
well-defined in the limit f ^ 0. 

Note that along a closed orbit, the phase difference of 
two interfering trajectories is fixed. 
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but B{t,t) = 2r > 0. 

Alternatively, one might use the implicit pair of equations 
l|21|l and l|28|l for this purpose. 



